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Abstract 

We use a statistical-thermodynamic model to investigate the formation and 
composition of noble-gas-rich clathrates on Pluto's surface. By considering 
an atmospheric composition close to that of today's Pluto and a broad range 
of surface pressures, we find that Ar, Kr and Xe can be efficiently trapped in 
clathrates if they formed at the surface, in a way similar to what has been 
proposed for Titan. The formation on Pluto of clathrates rich in noble gases 
could then induce a strong decrease in their atmospheric abundances relative 
to their initial values. A clathrate thickness of order of a few centimeters 
globally averaged on the planet is enough to trap all Ar, Kr and Xe if these 
noble gases were in protosolar proportions in Pluto's early atmosphere. Be- 
cause atmospheric escape over an extended period of time (millions of years) 
should lead to a noble gas abundance that either remains constant or in- 
creases with time, we find that a potential depletion of Ar, Kr and Xe in the 
atmosphere would best be explained by their trapping in clathrates. A key 
observational test is the measurement of Ar since the Alice UV spectrometer 
aboard the New Horizons spacecraft will be sensitive enough to detect its 

abundance ~10 times smaller than in the case considered here. 
Keywords: Pluto, surface - Pluto, atmosphere - Ices - Triton - 

Trans-neptunian objects 
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i 1. Introduction 



2 Gas hydrates, or clathrates, may exist throughout the solar system. Com- 

3 parison of predicted stability fields of clathrates with conditions in various 

4 planetary environments suggest that these structures could be present in the 

5 Martian permafrost (Musselwhite & Lunine 1995; Thomas et al. 2009; Swin- 
e die et al. 2009; Herri & Chassefiere 2012; Mousis et al. 2013), on the surface 
7 and in the interior of Titan (Tobie et al. 2006; Mousis & Schmitt 2008), and 
s in other icy satellites (Prieto-Ballesteros et al. 2005; Hand et al. 2006). It 

9 has also been suggested that the activity of many comets could result from 

10 the destabilization of these ices (Marboeuf et al. 2010, 2011, 2012a). On 

11 Earth, the destabilization of significant masses of CO2 and CH4 potentially 

12 stored in clathrates buried in seabeds and permafrost is regarded as a pos- 

13 sible aggravating factor in future global warming (clathrate gun hypothesis 

14 - Kennett et al. 2003). Broadly speaking, clathrates are thought to have 

15 taken part in the assemblage of the building blocks of many bodies of the 

16 solar system and may be in other planetary systems (Lunine & Stevenson 

17 1985; Mousis et al. 2002, 2006, 2009, 2010, 2011, 2012a; Mousis k Gautier 
is 2004; Alibert et al. 2005; Marboeuf et al. 2008; Madhusudhan et al. 2011; 

19 Johnson et al. 2012). 

20 Clathrates have also been proposed to be at the origin of the noble gas 

21 deficiency measured in situ by the Huygens probe in the atmosphere of Ti- 

22 tan (Osegovic & Max 2005; Thomas et al. 2007, 2008; Mousis et al. 2011). 

23 In the case of Mars, important quantities of argon, krypton and xenon are 

24 believed to be trapped in clathrates located in the near subsurface and their 

25 storage in these structures could explain the measured two order of magni- 
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26 tude drop between the noble gas atmospheric abundances in Earth and Mars 

27 (Mousis et al. 2012b). Here we investigate the possibility of formation of 

28 clathrates rich in noble gases on Pluto's surface. To do so, we use the same 

29 statistical-thermodynamic model applied to the case of Titan to determine 

30 the composition of clathrates that might form on Pluto and to investigate 

31 the possible consequences of their presence on the atmospheric composition. 

32 2. The statistical thermodynamic model 

33 To calculate the relative abundances of guest species incorporated in a 

34 multiple guest clathrate (hereafter MG clathrate) at given temperature and 

35 pressure, we use a model applying classical statistical mechanics that relates 

36 the macroscopic thermodynamic properties of clathrates to the molecular 

37 structure and interaction energies (van der Waals & Platteuw 1959; Lu- 

38 nine & Stevenson 1985). It is based on the original ideas of van der Waals 

39 and Platteeuw for clathrate formation, which assume that trapping of guest 

40 molecules into cages corresponds to the three-dimensional generalization of 

41 ideal localized adsorption. 

42 In this formalism, the fractional occupancy of a guest molecule K for a 

43 given type t (t — small or large) of cage (see Sloan 1998; Sloan & Kohn 2008) 

44 can be written as 

45 where the sum in the denominator includes all the species which are present 

46 in the initial gas phase. Cx,t is the Langmuir constant of species K in the 

47 cage of type t, and Pk is the partial pressure of species K. This partial 
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48 pressure is given by Pk = %k x P (we assume that the sample behaves 

49 as an ideal gas), with xk the mole fraction of species K in the initial gas, 

50 and P the total atmospheric gas pressure, which is dominated by N 2 . The 

51 Langmuir constant depends on the strength of the interaction between each 

52 guest species and each type of cage, and can be determined by integrating 

53 the molecular potential within the cavity as 

54 where R c represents the radius of the cavity assumed to be spherical, ks the 

55 Boltzmann constant, and WK,t{ r ) is the spherically averaged Kihara potential 

56 representing the interactions between the guest molecules K and the H 2 

57 molecules forming the surrounding cage t. This potential w(r) can be written 

58 for a spherical guest molecule, as (McKoy & Sinanoglu 1963) 
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(4) 



eo In Eq. 3, z is the coordination number of the cell. This parameter depends 

61 on the structure of the clathrate (I or II; see Sloan & Koh 2008) and on the 

62 type of the cage (small or large). The Kihara parameters a, a and e for the 



63 molecule- water interactions, given in Table A.l, have been taken from the 
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64 recent compilation of Sloan & Koh (2008) when available and from Parrish 

65 & Prausnitz (1972) for the remaining species. 

ee Finally, the mole fraction fx of a guest molecule K in a clathrate can be 
67 calculated with respect to the whole set of species considered in the system 
ea as 

, b s y K)S + b e y K ,e ,r\ 

j r ~ 1 u ' 

b s2^jyj,s + be)^jyj,e 

69 where b s and bi are the number of small and large cages per unit cell respec- 

70 tively, for the clathrate structure under consideration, and with ^2 k Jk = 1- 

71 Values of R C) z, b s and bi are taken from Parrish & Prausnitz (1972). 

72 In the present approach, the dissociation pressure of the MG clathrate 

73 and the mole fractions of the trapped volatiles are independently calculated. 

74 All mole fraction calculations are performed at the dissociation pressure P = 

75 P^tx °f tli e clathrate, i.e. temperature and pressure conditions at which this 

76 ice forms. This dissociation pressure can be deduced from the dissociation 

77 pressure P* ss of a pure clathrate of species K, as (Hand et al. 2006; Thomas 

78 et al. 2007): 

-^mix = f pdiss) ' (6) 

K K 

79 where x K is the atmospheric mole fraction of species K and P^ ms its disso- 
so ciation pressure. P^ 1SS derives from laboratory measurements and follows an 
si Arrhenius law (Miller 1961) as 



log(P diss ) = A + |, (7) 
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82 where P dlss is expressed in Pa and T is the temperature in K. The constants 

83 A and B for N2 and Ar have been fitted to experimental data (Lunine & 

84 Stevenson 1985; Sloan 1998) and those for Xe and Kr have been taken from 



85 Fray et al. (2010) and Marboeuf et al. (2012b), respectively (see Table A. 2) 



3. Clathrate production on Pluto 



87 



Depending on its composition, the dissociation pressure of MG clathrate 

as varies between ~5.1 x 10" 10 and 2.1 x 10" 5 Pa at Pluto's average surface 

89 temperature (~50 K; Lellouch et al. 2000), a value well below the atmo- 

90 spheric surface pressure which lies in the 0.65-2.4 Pa range (Elliot et al. 

91 2003; Sicardy et al. 2003). This implies that MG clathrate remains stable 

92 at the surface irrespective of Pluto's seasonal variations and that its dissoci- 

93 ation and reformation cannot occur under the planet's current atmospheric 

94 conditions. The only conditions on Pluto's surface allowing MG clathrate 

95 formation from an initial inventory present in the atmosphere require that 

96 the surface be hotter than the present-day temperature, most likely at early 

97 epochs after the planet's formation (see e.g. Fig. 1 of McKinnon 2002) or 

98 during the collisions that engendered the Pluto-Charon binary system (McK- 

99 innon 1988, 1989). An alternative possibility would be the recent or ancient 

100 release of hot ice from the interior of Pluto as the result of cryovolcanic 

101 events (Cook et al. 2007). It is during clathrate formation that substantial 

102 amounts of volatiles might have been sequestrated from the atmosphere into 

103 the surface. 



Table [AT3] gives the composition of Pluto's proto-atmosphere used in our 



105 calculations. We made the conservative assumption that all noble gases were 
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106 initially present in the proto- atmosphere of Pluto, with Ar/N, Kr/N and 

107 Xe/N ratios assigned to be protosolar (Asplund et al. 2009). Because Pluto's 
los proto-atmosphere is expected to be strongly dominated by N 2 and that N 2 
109 clathrate is of structure II (Lunine & Stevenson 1985; Sloan & Koh 2008), 
no we show here calculations of the MG clathrate composition only for this 
in structure. It is also important to note that our composition calculations are 
n2 only valid along the dissociation curve of the clathrate of interest (see Fig. 



A.l). 



H4 From our calculations, we find that three kinds of clathrates with distinct 
us compositions might form on Pluto's surface, each of them containing noble 



lie gases in different proportions. Figure A. 2 shows the composition of these 

in clathrates computed for an atmospheric pressure ranging between 1 and 10 3 

us Pa. Note that these surface pressures correspond to a MG clathrate equi- 

n9 librium temperature in the ~80-100 K range. If a greater surface pressure 

120 is considered, then the MG clathrate will also form at a higher equilibrium 

121 temperature. A first layer forms from the gas phase composition depicted in 

122 Table |AT3 Irrespective of the pressure considered, the mole fraction of Xe 

123 trapped in this clathrate is between ~0.15 and 0.76, i.e. a range of values 

124 that is ~31, 000-159, 000 times larger than its atmospheric mole fraction. The 

125 mole fraction of Kr is also substantially enhanced by a factor of ~400-750 in 

126 clathrate compared to its atmospheric value. In contrast, the mole fraction 

127 of Ar trapped in clathrate evolves from a slight impoverishment (~0.5) to a 

128 moderate enrichment (~3) with increasing pressure, compared to its atmo- 

129 spheric value. The second clathrate layer forms once Xe is fully trapped in 

130 the first layer. In this case, only N 2 , Kr and Ar remain in the gas phase (the 
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131 mole fractions of these species trapped in the first layer remain negligible). 

132 The mole fraction of Kr trapped in this clathrate is enhanced by a factor 

133 of 400-4,000 times compared to its atmospheric mole fraction in the 1-10 3 

134 Pa pressure range. The Ar mole fraction in this layer is also found to be 

135 moderately enriched by a factor of ~2-3, compared to its atmospheric mole 

136 fraction. In its turn, a third clathrate layer forms when Kr is fully trapped 

137 in the second layer and in this case only N 2 and Ar remain in the coexisting 

138 gas phase. The fraction of Ar trapped in this clathrate remains constant 

139 irrespective of the surface pressure considered and is found to be ~3 times 
wo larger than its atmospheric value. 

141 If Ar, Kr and Xe were initially in protosolar abundances in the atmosphere 

142 of Pluto, the amount of clathrates needed for their sequestration is relatively 

143 low. For example, if the three clathrate layers formed at their equilibrium 

144 temperature^] for a surface pressure of 2.4 Pa, their total equivalent thickness 

145 is of order ~4 mm globally averaged on the planet, assuming a full clathration 

146 efficiency and the presence of a structure II clathrate. Interestingly, calcula- 

147 tions conducted in the case of formation of structure I clathrate lead to similar 

148 conclusions. The noble gas trapping efficiencies are still very high but require 

149 an overall thickness of the three clathrate layers of ~2 cm. In both clathrate 

150 structures, the equivalent layer of the Ar-dominated clathrate is more than 3 

151 orders of magnitude thicker than those of Xe- and Kr-dominated clathrates. 

152 One must also note that the noble-gas-rich clathrates formed on Pluto could 

153 consist in a mixture of structures I and II clathrates. Indeed, in our com- 

lr The equilibrium temperatures of the first, second and third clathrate layers are ~77, 
76 and 76 K, respectively. 
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154 
155 
156 



putations of the composition of structure II clathrates, we find that the first 
layer formed is dominated by Xe, which is itself predicted to form a structure 
I clathrate (Sloan & Koh 2008). 



is? 4. Competition with atmospheric escape 

158 An alternative method for losing these noble gases from Pluto's atmo- 

159 sphere could be atmospheric escape. The atmosphere of Pluto is expected to 
wo escape efficiently due to Pluto's low gravity (Strobel 2008), though debate 
lei exists as to whether the escape rate is greater than subsonic (Tucker et al. 

162 2012). If the escape of Pluto's atmosphere is hydrodynamic (Trafton et al. 

163 1997; Strobel 2008), the outflow of N 2 provides enough energy to drag Ar, 

164 Kr and Xe from the atmosphere. The escape rate of these gases depends 

165 directly on the escape rate of N 2 (Hunten et al. 1987): 

^ X; „ / 17l c — m* \ 

F t = -^F N2 — ? l - , 8 

X N2 \m c -m N2 J 

lee where F is the escape flux, X is the abundance of the noble gas (subscript i) 

167 and N 2 , m is the mass of each constituent and m c is the critical mass which 

168 is also a function of the N 2 escape rate (Hunten et al. 1987): 

, kTF N2 

m c = m N2 + -— , (9) 

b g X N2 

leg where k is the Boltzmann's constant, T is the temperature, b is the binary 

170 diffusion coefficient given as b = AT S with A and s constants determined 

171 through modeling or laboratory measurements, and g is the acceleration due 

172 to gravity. 
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173 It is clear from Eq. [9] that the critical mass is always greater than the 

174 mass of N2, so constituents with mass less than N2 are always subject to 

175 escape when hydrodynamic escape is occurring. Ar, Kr and Xe all have 

176 masses greater than N 2 , but the proposed escape rate of N 2 is so high that 

177 the critical mass is several times greater than the mass of Xe and the last term 

178 in Eq. [8] has a value of 1.0. This means that the escape rate of the noble gases 

179 relative to the escape rate of N 2 , under hydrodynamic escape conditions, is 
wo equivalent to the abundance of the noble gas relative to N 2 as suggested in 

181 Eq. [8j Therefore, at the highest possible escape rate for the noble gases, 

182 the abundance of the noble gas relative to N 2 will remain constant. In the 

183 case of escape rates that are subsonic, escape of the noble gases will be far 

184 less efficient due to insufficient energy to drag the heavier molecules from 

185 the atmosphere. As a result, the noble gas escape rate relative to N 2 's rate 
we of escape will be less than the abundance of the noble gases relative to N 2 
187 in the atmosphere. Over time, the differential escape rates will lead to an 
we increase in the noble gas abundance in Pluto's atmosphere. This means that 

189 without trapping of Ar by clathrates, escape over an extended period of time 

190 (millions of years) will lead to a noble gas abundance that either remains 

191 constant or increases with time. Therefore, a depletion of Ar, Kr and Xe 

192 in the atmosphere can best be explained by trapping of the noble gases in 

193 clathrates. 



194 5. Discussion 



195 Using the noble gas abundances given in Table A.3 and the fitting laws 
we of sublimation laboratory data proposed by Fray & Schmitt (2009), we find 
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197 that the equilibrium temperatures of Ar, Kr and Xe pure ices are ~39.5, 

we 42.5 and 55.3 K at the surface pressure of 2.4 Pa, respectively. Given Pluto's 

199 average surface temperature (~50 K), Xe appears stable on the surface in 

200 present-day conditions and the condensation of this noble gas should also 

201 induce a decrease of its atmospheric abundance. This implies that there is 

202 no way to disentangle Xe's previous trapping in clathrate versus its simple 

203 condensation. This is fortunately not the case for Ar and Kr, so the mea- 

204 surement of at least one of these two noble gases is critical for testing the 

205 scenario of clathrate trapping. The UV spectrometer Alice aboard the New 

206 Horizons spacecraft will be sensitive enough to detect an argon abundance 

207 that is ~ 10 times smaller than protosolar Ar/N in some circumstances (see 



208 appendix Appendix A). 

209 Our scenario of the noble gas trapping in clathrates on Pluto is motivated 

210 by the fact that the same interpretation was provided in the case of Titan 

211 in order to account for its observed noble gas deficiency (Osegovic & Max 

212 2005; Thomas et al. 2007, 2008; Mousis et al. 2011). The situation here 

213 is much more favorable than in the case of Titan since the clathrate layer 

214 required on Pluto's surface is ~10 4 times thinner. This is essentially due to 

215 the difference between the surface pressures since the two atmospheres are 

216 both dominated by N2. 

217 Note that our model considers only the trapping of the noble gases that 

218 were present in the early atmosphere. It does not exclude the possibility 

219 that the bulk of these noble gases is still in the interior of Pluto. However, 

220 if there is outgassing, the noble gases released in the atmosphere should be 

221 trapped as well by clathrates. Our calculations are also valid irrespective of 
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222 the source (primordial or radiogenic) of the noble gases potentially present 

223 in Pluto's atmosphere. Indeed, since Pluto is ~half rock, a significant part 

224 of the existing argon could result from the radiogenic decay of potassium-40. 

225 Interestingly, similar calculations have been performed in the case of Tri- 

226 ton and they lead to conclusions as favorable as in the case of Pluto. However, 

227 the surface temperature of Triton, ~38 K (Tryka et al. 1994), is much lower 

228 than Pluto's. It is low enough that the condensation of the three noble gases 

229 as pure ices on the surface would largely remove them from the atmosphere 

230 irrespective of their sequestration in clathrate. Nonetheless, a mass spec- 

231 trometer directly sampling Triton's atmosphere could in principle detect the 

232 atmospheric abundances of all three noble gases to assess if they were con- 

233 sistent with coexistent surface ices. Since the atmospheric abundances in 

234 coexistence with crustal clathrate are orders of magnitude smaller (probably 

235 too small for detection), this provides an eventual test of the hypothesis in 

236 the case of a future mission to the Neptune system. 

237 6. Conclusions 

238 By considering an atmospheric composition close to that of today's Pluto 

239 and a broad range of surface pressures, we find that Ar, Kr and Xe can be 

240 efficiently trapped in clathrates if they formed at the surface. The formation 

241 of noble gas-rich clathrates on Pluto could then induce a strong decrease of 

242 their initial atmospheric abundances. A clathrate thickness of order of a few 

243 centimeters globally averaged on the planet is indeed enough to trap Ar, Kr 

244 and Xe irrespective of the clathrate structure, if they were in protosolar pro- 

245 portions in the early Pluto's atmosphere. We suggest that the measurement 
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246 of the Ar abundance by the Alice ultraviolet spectrometer (Stern et al. 2008) 

247 aboard the New Horizons spacecraft during Pluto's flyby should provide a 

248 test of the validity of our scenario. 

249 Acknowledgements 

250 O. Mousis acknowledges support from CNES. J.I. Lunine was supported 

251 by a contract from JPL under the Distinguished Visiting Scientist program. 

252 J. H. Waite acknowledges support from NASA Rosetta funding. 

253 Appendix A. Estimating Argon Airglow Emission for the New 

254 Horizons Mission 

255 With the nominal Pluto encounter of the New Horizons mission fully 

256 planned, we wish to estimate the amount of Ar I that can be detected by the 

257 Alice UVS instrument during approach. There will be 11 Alice observations 

258 at close range on approach dedicated to measuring extended airglow emission, 

259 with integrations times ranging from 360 to 6900 seconds each, for a total of 

260 approximately 6 hours of integration. The midpoint distances from Pluto for 

261 these observations range from 1 x 10 6 km to 3 x 10 5 km. We use the planned 

262 observation times and their corresponding midpoint distances, as well as the 

263 Alice effective area at the Ar I doublet (104.8, 106.7 nm), to estimate the 

264 signal to noise ratio (SNR) for a given Ar I mixing ratio. 

265 The Ar I brightness was determined using the Atmospheric Ultraviolet 

266 Radiance Integrated Code (hereafter AURIC) (Strickland et al. 1999). This 

267 was run with Ar I abundances of 0.1%, 0.3%, 1%, 3%, and 10%, to pro- 

268 ducc volume production rates from photoelectron impact and photoexcita- 
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269 tion processes. "Model 2" densities from Krasnopolsky & Cruikshank (1999) 

270 were used for N2 and CH4 in the model atmosphere (7.1 x 10 13 cm -3 and 

271 6.4 x 10 11 cm -3 , respectively), with laboratory measured cross sections and 

272 TIMED/SEE solar spectral irradiance for excitation. The Ar I brightnesses 

273 fed into the Alice UVS model thus are the emergent intensities upwelling 

274 through the N2/CH4 atmosphere to space. 

275 Figure A.3| shows simulated brightness images for one of the approach 



276 observations, along with the observation coordinates. Figure A. 4 shows the 



277 curve of growth made with the brightness estimates for each Ar I abundance 

278 simulated. We use this brightness relation to calculate the SNR obtained with 

279 the Alice airglow observations. The total counts obtained (C in photons) in 

280 the Alice slit is just: 

10 6 

C = — R A (A.l) 

281 where R is the brightness (Rayleighs), and A is the effective area (0.1 cm 2 at 

282 the Ar I doublet). Since each observation differs in length (ti in seconds) and 

283 is performed at different distances (leading to different solid angles), we sum 

284 the photons piecewise. f2, is the effective solid angle of the emitting object 

285 visible in the slit. The different exposure times (ti), midpoint distances from 

286 Pluto (di), and slit-filling solid angle for each observation are listed in Table 

287 I A. 41 

288 For our SNR estimates we consider only the emission from the area of 

289 Pluto's solid disk. We assume a detector background dark rate of 0.02 counts 

290 per spatial/spectral element (as measured in flight) and the AURIC bright- 

291 ness estimates of the 104.8 and 106.7 nm Ar I emission lines to calculate the 
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292 SNR for a given mixing ratio, shown in Figure A. 5 



293 With these calculations we conclude that in 6 hours of observing on final 

294 approach to Pluto, the Alice instrument should be able to detect Argon emis- 

295 sion in Pluto's atmosphere even at lower than solar abundances. Assuming a 

296 number abundance of 0.25 ppm for Ar I and 36 ppm for N2, i.e., at one tenth 

297 solar abundance, the SNR will be ~16. Unfortunately, the measurements of 

298 Kr and Xe abundances in the atmosphere of Pluto is beyond the instrument's 

299 capabilities. 
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Table A.l: Parameters for Kihara and Lennard- Jones potentials 



Molecule ctr-w (A) e k-w/^b (K) a^_iy (A) Reference 



N 2 


3.0993 


133.13 


0.3526 


Herri & Chassefiere (2012) 


Ar 


2.9434 


174.14 


0.184 


Herri & Chassefiere (2012) 


Kr 


2.9739 


198.34 


0.230 


Parrish & Prausnitz (1972) 


Xe 


3.32968 


193.708 


0.2357 


Sloan & Koh (2008) 



&k-w is the Lennard-Jones diameter, 6k-w is the depth of the potential 
well, and ax~w is the radius of the impenetrable core, for the guest-water 

pairs. 
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Table A. 2: Parameters of the dissociation curves for various single guest clathrate hydrates. 
A is dimensionless and B is in K. Constants for Kr and Xe are given for an Arrhenius law 
making use of Napierian logarithm. 



Molecule 


A 


B 


N 2 


9.86 


-728.58 


Ar 


9.34 


-648.79 


Kr 


22.3934 


-2237.82 


Xe 


16.62 


-3159 
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Tabic A. 3: Assumed composition of Pluto's atmosphere at the ground level. Noble gas 
abundances relative to N 2 are assumed to be protosolar (Asplund et al. 2009). 



Species K Mole fraction %k 



N 2 


9.31 


X 


io- 1 


Ar 


6.92 


X 


io- 2 


Kr 


4.90 


X 


10~ 5 


Xe 


4.79 


X 


10~ 6 
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Tabic A. 4: Relevant parameters of each airglow observation. 



Observation 


di (km) 


U (s) 


Qi (/xrad) 


PC-Airglow-Fill-0 


1116455 


6900 


3.21 


PC-Airglow-Fill-2 


819350 


600 


4.70 


PC-Airglow-Fill-2 


809012 


600 


4.77 


PC-Airglow-Fill-2 


798261 


600 


4.84 


PC-Airglow-Fill-2 


788336 


600 


4.91 


PC-Airglow-Fill-2 


780204 


360 


4.97 


P C-Airglow- Appr- 1 a 


611833 


1200 


6.49 


PC-Airglow-Appr-lb 


575456 


600 


6.94 


P C-Airglow- Appr-2 


525813 


4200 


7.63 


P C-Airglow- Appr-3 


408242 


3600 


9.93 


P C-Airglow- Appr-4 


318118 


3000 


12.83 
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Pressure (Pa) 

Figure A.l: From top to bottom: equilibrium curves of Xe-, Kr-, and Ar-dominated 
clathrates (the latter two appear superimposed) . The arrow pointing down represents the 
path followed by a cooling atmosphere with a surface pressure of 2.4 Pa. When the cooling 
curve intercepts an equilibrium curve, then the corresponding clathrate forms. 
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Figure A. 2: From top to bottom: mole fractions of volatiles encaged in clathrates that 
successively form and calculated as a function of the surface pressure of N2. Noble gas 
abundances are assumed to be protosolar relative to N 2 (Asplund et al. 2009). The 
clathrate composition is investigated from a gas phase composed of Ar, Kr, Xe and N2 
(top panel), of Ar, Kr and N2 (middle panel), and of Ar and N2 (bottom panel). 
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Figure A. 4: Emergent Ar I 104.8 and 106.7 nm brightness as a function of mixing ratio. 
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Figure A. 5: Signal to noise ratio for the added Ar I 104.8 and 106.7 nm line emission vs. 
mixing ratio. 
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